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Abstract. I review recent progress in analysing deep inelastic scattering structure 
functions in global analyses. The new ingredients are new data and attempts to incor¬ 
porate heavy quarks consistently. A new way of including the resummation of large 
\ogl/x terms is discussed. 


INTRODUCTION 

Since we last met in Rome last year there has been progress in onr understanding 
of deep inelastic scattering (DIS) structure functions and the implications of new 
measurements in the context of global analyses and extraction of parton distribu¬ 
tion functions (pdf’s). Firstly there have been updates to two major experiments 
— the muon-nucleon NMC collaboration at CERN and the neutrino-iron CCFR 
collaboration at Fermilab. The NMC [1] has produced their final numbers for F 2 
for proton and deuterium targets at 90, 120, 200 and 280 GeV while the CCFR [2] 
collaboration has implemented new energy calibrations to its analysis. The latter 
data are consistent with a greater value of a s (M z ) than the previously reported 
value. Finally there have been the first reported measurements from HERA on the 
charm structure function Ff [3,4] which indicate that a substantial fraction of the 
total F 2 at HERA comes from charm production. 

On the theoretical side, progress has been made following improved treatments 
of heavy quark production in DIS. Why is this important? From above we see 
that the component of F 2 is relatively large and so it is obviously important to 
have a consistent description of this component and to check this with boson-gluon 
fusion mechanisms. Also if we are trying to understand the small x behaviour of F 2 
(i.e. BFKL versus GLAP) then the data force us to study small Q 2 i.e. the region 
where we move through the charm threshold; so this should be understood. Finally 
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it is necessary to produce pdf’s for charm and bottom flavours to insert into other 
processes such as jet production. 

Large x. high Q 2 is where the excitement resides just now. Can we be precise 
about the conventional background and estimate the uncertainty in the cross section 
around Q 2 ~ 10 5 GeV 2 , x ~ 0.45? I shall discuss this briefly. 

For some time it has been realised that in addition to resumming large logs in Q 2 
(GLAP), there may be significant contributions arising from resumming potentially 
large logs in 1/x (BFKL). The difficulty is to devise a theoretically consistent 
procedure of ‘ marrying’ the two in a practical way that confronts the data and 
tries to answer the question whether these data favour or disfavour the inclusion 
of the log 1/x resummed terms. Here there has been definite progress and I shall 
highlight a recent analysis by Thorne [5]. 


RECENT GLOBAL ANALYSES 

The two major providers of pdf’s, CTEQ and MRS, continue to update their 
global analyses of DIS and other processes. Both are now attempting to include a 
theoretically improved treatment of charm and bottom and these treatments are 
discussed in the following section. In this section I will highlight some results 
emerging from these on-going analyses. Lai and Tung [6] demonstrate the effect 



FIGURE 1. Comparison of HI data [8] with calculations using CTEQ4M pdf’s in the original 
massless scheme (solid line) with the new fit CTEQ4HQ using the ACOT procedure (dotted line). 

of including the variable flavour scheme (VFNS) developed by Aivazis et al. [7] 
(ACOT) based on their ealier work [9] into the CTEQ analysis. This is to be 





contrasted with the previous treatment in CTEQ where charm evolved as a massless 
quark. The effect is illustrated in fig. 1. 

Only a small change is observed as one changes from CTEQ4M (the massless 
evolution) to the improved VFNS procedure which is more noticeable at small x. 
A comparison of the resulting pdf’s at Q 2 = 25 GeV 2 is shown in fig. 2 Moreover, 



FIGURE 2. Comparison of CTEQ4M and CTEQ4HQ pdf’s 

the quality of the resulting fit to the full F 2 is improved at low x. 

The series of global analyses carried out by MRS has been updated. The new 
data from NMC [1] and CCFR [2] are included and the improved treatment of the 
heavy quark developed by Martin et al. [10] is used (see next section). As with 
the previous MRS analysis [11] a low value of the starting scale Ql of 1 GeV 2 is 
taken. This is mainly to reach as small x as possible for the HERA data [8,12] 
for which a cut at Q 2 = 1.5 GeV 2 is taken while the cut is at 2 GeV 2 for other 
DIS data. In order to study the sensitivity to the value of a s (M z ), the analysis is 
repeated at fixed values of a s (M z ) over a wide range. The resulting values of the 
individual y 2 values for the DIS data sets are shown in Eg. 3. From the figure we 
see (i) that the HERA [8,12] data prefer a value of a s (M z ) around 0.12 or so [13], 
(ii) that the updated NMC [1] data also prefer this somewhat larger value and (iii) 
the re-analysed CCFR [2] F 2 data strongly favour a value at the upper end. This 
leaves BCDMS [14] data as the one set which still favours a relatively low value of 
a s (M z ). But we should recall that these data are very precise and measure scaling 
violations in the large x region where the determination of a s (M z ) is insensitive to 
possible uncertainty of the gluon distribution. 

The resulting best value of a s (M z ) which emerges from this overall global analysis 
is 0.118 and we show comparisons of the fit using this value with a selection of data. 
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FIGURE 3. Values of y 2 for various DIS datasets [8,12,14,1,15,16] as a s (Mz) is varied 

The small x data comparison is shown in fig. 4. The comparsion with the new NMC 
data is shown in fig. 5. Notice that points below Q 2 = 2 GeV 2 , while shown in the 
figure, were not included in the fit. Finally we show the comparison with the new 
data from CCFR [2] on F 2 extracted from v and v interactions off an iron target 
in Fig. 6. The failure to agree below x of around 0.1 is well known. Basically 
these low x data prefer not to have the shadowing corrections, observed in muon 
DIS, applied. Because of this conflict (with the NMC data for example) CCFR F 2 
below x = 0.1 is dropped in the fit. Comparing the curves and large x data in fig. 
6 one can see that the data prefer slightly stronger Q 2 variation than that from 
using a s (M z ) = 0.118. Much better consistency is found for xF 3 across the entire 
Q 2 , x range. 

In this analysis, the pdf’s at Ql are parametrised in the usual MRS [11] way 
with xg(x, Qq) ~ x~ Xg and xS(x,Qq) ~ x~ Xs for the glue and sea as x — > 0. The 








FIGURE 4. Comparison of the new fit using a s (Mz) = 0.118 with the small x data from refs 
[8,12,16]. 


exponents X g and X s are not constrained and as in the previous analysis one finds 
that at such a low choice of Qq — 1 GeV 2 X s is positive (singular) while X g comes 
out negative — i.e ‘valence-like’ gluon. The actual value of the exponent varies for 
different values of a s (Mz) (—0.15 to —0.21) but when Q 2 reaches just 2 GeV 2 all 
the gluon pdf’s are approximately ‘flat’, see fig. 7. Taking the fit with a s = 0.118 
and evolving to Q 2 = 25 GeV 2 yields pdf’s which are extremely close to those 
shown in fig. 2 so the different treatments of charm partons in the latest analyses 
of CTEQ and MRS actually lead to little numerical difference in the resulting pdf’s. 





FIGURE 5. Comparison of the new fit using a s (Mz) = 0.118 with the NMC data [1]. 


THEORETICAL DISCUSSION OF HEAVY QUARK 
PRODUCTION IN DIS 

Let us consider possible approaches to treating the charm production contribu¬ 
tion to the structure function F 2 . 

Massless parton evolution 

The most simplistic approach is to assume that a probe of virtuality Q 2 can 
resolve a charm quark pair in the proton sea when Q 2 > m 2 . Since such pairs 
originate from fluctuations of the gluon field, g —> cc, a perturbative treatment 
should be valid as long as m 2 3> A q CD - As Q 2 increases, 0(m 2 /Q 2 ) corrections 
to the standard GLAP evolution become less important, and the charm quark can 
be treated as a (fourth) massless quark. These ideas are embodied in the ‘massless 
parton evolution’ (MPE) approach 


c(x, Q 2 ) = 0 for Q 2 < n 2 c , 






FIGURE 6. Comparison of the new fit using a s (Mz) = 0.118 with the new CCFR [2] data on 
F‘ 2 . Here the theory curves have been multiplied by a heavy target x dependent correction factor 
and only the statistical errors are shown on the data. 


Uf = 3 + 6(Q' 2 - l4) in Pqg, Pgg, ft),... , (1) 

where /i c = 0(m c ). The charm contribution to the structure function is then 

F£( x , Q 2 ) = | xc(x, Q 2 ) (2) 

in lowest order. This is the approach adopted at NLO in the MRS and CTEQ 
global parton analyses, with /i c chosen in the MRS analysis to achieve a satisfactory 
description of the EMC F| data [17], while in the CTEQ analysis /i c is set equal 
to m c . For example, in the MRS(A) analysis [18] it was found that n 2 = 2.7 GeV 2 
and that this was to a good approximation equivalent to taking 

2c(x, Ql) = SS(x,Q 2 0 ) 


( 3 ) 





FIGURE 7. Gluon distribution functions at Q 2 = 2 GeV 2 as a s (Mz) is varied from 0.105 
(upper dashed) through 0.110 (dot-dashed), 0.115 (upper dotted), 0.1175 (upper solid), 0.120 
(lower dashed), 0.125 (lower solid) and 0.1275 (lower dotted). 


with 5 « 0.02 at the input scale Qq = 4 GeV 2 . That is at the input scale, charm 
(c + c) was found to have approximately the same shape as the total quark sea 
distribution S , and moreover to form about 2% of its magnitude. 

In NLO the charm structure function is given by 


f;(x,q 2 ) = l I 1 dz 
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while for the gluon we have 
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in (4) has the form 
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and this expression with the standard massless MS coefficient functions was used 
to calculate F!> in the MPE approach. 




Boson-gluon fusion 


A procedure which takes a very different approach is to generate F 2 C from photon- 
gluon fusion (PGF) [19,20] alone. In this case, the charm pdf appearing in eq.(4) is 
identically zero everywhere but the gluon coefficient function involves the non-zero 
charm mass m c and is the well-known PGF cross section, i.e. C g (z) —> C g ¥ (z, Q 2 ) 
with 
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where f3 is the velocity of one of the charm quarks in the photon-gluon centre-of- 
mass frame 
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( 8 ) 


The 0 function in eq.(7), Q(W 2 — 4 ml), represents the cc production threshold, 

where W is the c.m. energy. Its presence guarantees (3 2 > 0. The exact next-to- 

leading order corrections to the PGF structure function are known [21] and recently 

[ 22 ] an analysis has been used to perform [as ln(( 3 2 /m[[)) n resummation up to 
2 

0[^j) for Q 2 S> m 2 c . Such an approach is of course not applicable in the threshold 
region Q 2 > m 2 . However for many purposes the PGF procedure is insufficient. If 
we consider the high Q 2 limit of we get 


C 'J GF ^ {z 2 + (1 - z) 2 } In + 82(1 -z)-l (9) 


as m c —> 0 , which differs from the exact m c = 0 coefficient by the presence 
of Q 2 /m 2 in the argument of the logarithm. This is a signal of potentially large 
logs which should be resummed by GLAP evolution. Thus at large Q 2 we should 
include the charm quark as a parton in GLAP evolution. In fact we should be able 
to provide a charm pdf all the way down to the charm threshold in order to insert 
all relevant pdf’s into calculations of other perturbative QCD processes. The goal 
therefore is to do this in a consistent way 


Variable flavour number scheme procedure 

Here the procedure is to have no charm pdf below Q 2 = m 2 but for Q 2 > m 2 to 
have c + c ~ a s P qg In Q 2 /m 2 , where P qg is the splitting function g —> cc evaluated 
with rn — 0. Below Q 2 = m 2 charm is produced solely by boson-gluon fusion 



as above and this is labelled the ‘hxed flavour number’ scheme (FFNS). Thus for 
rif = 3 we have an OPE for the structure functions 

F = C F fns ~jjt ’ w ( h 2 )^ 0 FFN s(a s (fi 2 )) + O(-) ( 10 ) 

Above Q 2 = rri 2 we work in the ‘variable flavour number’ scheme (VFNS) with 
rif — 4 and write the OPE 

fn 2 \ ~ O 2 m 2 in 2 

F = C VFNS Oi s ([i 2 ) O vfns (, a s ([i 2 )) + 0(- f) (11) 

\ h / /i /i V 

While the operators Ovfns nominally depend on m 2 they can be chosen so that 
they evolve according to m = 0 RGE’s. Thus in the VFNS the m 2 dependence ap¬ 
pears only in the corrections. The coefficients Cvfns are simply the usual massless 

MS coefficients C^g. The two schemes must be mutually consistent at all Q 2 > m 2 

2 

which implies (from eqns.(10,ll)) that the f) contribution in eq.(ll) can be 
written in the form 

rtr 2 ^ A 2 

C{-1) 0 FFNS (a s (fi 2 )) + O(-) (12) 

2 

where C' is a coefficient which ~ powers of Fjt. If we express the relation between 
the operators of the two schemes through a matrix A, 

Ovfns = A ^—-^,a s (fi 2 )^ Offns(c^s) (13) 

then, in particular, the charm pdf is determined purely in terms of the light parton 
pdf’s — principally the gluon distribution of course. The structure functions above 
Q 2 = m 2 can again be expressed in terms of the VFNS operators but now up to 
O(F^) corrections, 

(O 2 \ in 2 1 - A 2 

F= Cvfns f —, ot s j + C'{—f)A~ 2 Ovfns + C *(—) 

\ R / h \ O 

(O 2 m 2 \ ~ A 2 

= Cacot ( —, — f, 0! s ) Ovfns + £*(—) (14) 

\ [J, fj, J v 

where Cacot is the coefficient in the ACOT [7,9] formalism which then depends 
on m 2 . Thus in the ACOT formalism, while the partons obey a simple evolution 
(independent of m 2 ), the coefficients would have a complicated dependence on m 2 . 
From eqns (10,14) we obtain the relation between the coefficient functions 2 

2 i I am grateful to Robert Thorne for clarifying my understanding of the connection between the 
FFNS and VFNS. 
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where the term C'A^ 1 on the rhs —► 0 as Q 2 —> 00 . 

Buza ef al. [22] have evaluated the coefficient functions for Q 2 3> m 2 in the 
VFNS and the FFNS to NLO. They find that for Q 2 > 20GeV 2 the two schemes 
agree very closely. The same authors have worked out the relation (13) in the 
limit Q 2 3> ml exactly through terms ln 2 (m 2 //P), In(m 2 //P) and constants. This 
infers a matching of all the four flavour pdf’s above the scale g 2 = ml, to the three 
flavour pdf’s below. Again this is only for Q 2 3> m 2 , which does not solve the 
problem of what to do at Q 2 ~ m 2 c . So there has to be another matching between 
the calculated cc production cross section in the threshold region to the charm pdf 
prescription at Q 2 3> m 2 . In leading order this is exactly the ACOT procedure. It 
remains to be investigated how to do this retaining terms 0(m 2 c /Q 2 ) at NLO. 

Thus although the FFNS —> VFNS change of factorization scheme is very attrac¬ 
tive in principle, a major problem arises when one tries to go down in Q 2 (still at 
NLO) all the way to the charm threshold - retaining the C'A^ 1 term — which is the 
eventual aim of ACOT. The A -1 term is the one which contains the subtraction in 
the gluon coefficient functions for Fy in eq.(2) to avoid double counting which at 
LO is quite straightforward and takes the form 

cf\z) = Cg GF (z) - a 8 CW x pW In Q 2 /ml 

The effect of this is that P 2 C i s described at low Q 2 by PGF but as Q 2 —> 00 by the 
evolved charm pdf. 

Since the way it has been formulated is only at leading order (the charm pdf 
evolves only via P^\z)) this allows the standard expression to be used for the 
gluon coefficient function. We are still left with the problem of how to evolve 
all the pdf’s beyond L.O. - in the MS scheme say. While the splitting functions 
will continue to have the m — 0 forms, the entire m 2 dependence is put into the 
coefficient functions which will make them (the light parton coefficient functions 
as well as the charm coefficient functions) extremely complicated when they are 
eventually computed. 


MRRS procedure 

In view of the practical difficulties of implementing the ACOT procedure at NLO, 
Ryskin and MRS [10] introduce an alternative prescription which is consistent at all 
orders but which generates partons which may be used with the conventional MS 
coefficient functions. The charm quark evolution is modified however — of course 
only for Q 2 S> m 2 c does it evolve according to the massless MS. The modified split¬ 
ting function is derived by studying the leading log decomposition of the Feynman 
diagrams and takes the form 
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Thus as Q 2 3> m 2 , P cg (z, 73 #) —»> P qg (z) U=o in the MS scheme. The NLO 
splitting function obtained from the ladder graph actually has a natural scale given 
by Q 2 = (m 2 + k^)/z( 1 — z ). This choice was not the one chosen when the MS 
scheme was set up; instead the scale Q 2 = (m 2 + kj) was chosen and it is the latter 
scale which leads to the expression in eq.(16). The ‘natural’ choice of scale would 
have the advantage that it can be interpreted as corresponding to a ‘resolution 
threshold’ of Q 2 = Am 2 below which Q 2 is too small to allow sufficient time to 
observe the g —> cc fluctuations. To avoid double counting the gluon coefficient 
function for must have a subtraction term, i.e. 


where 
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(18) 


This subtraction ensures the required cancellation of the leading contribution 
to the charm quark evolution around Q 2 ~ m 2 , but the NLO contribution from 
the P£p term is still significant. However inserting the condition that the natural 
scale should really be Q 2 ~ 4m^ into the charm quark coefficient function forces 
the low Q 2 description of F 2 C to be given entirely by the PGF term. The resulting 
description of the measured cross sections is remarkably successful - see fig. 8 . 


LARGE X 

Because of the topical interest in the magnitude of the cross section at large 
x and large Q 2 it is worthwhile to study the precision by which we believe we 
understand the conventional DIS cross section. This translates into estimating our 
confidence in the large x pdf’s. In fig. 9 we show the electromagnetic contribution 
to F -2 at x = 0.45. For Q 2 < 200 GeV 2 the data have a normalisation uncertainty 
of about 3%. Taking the pdf’s from the latest global analyses together with a 
generous uncertainty in the value of a s (Mz) — from 0.110 to 0.125 we see that this 
gives an additional 7% to the extrapolated value at Q 2 ~ 10 5 GeV 2 . Including the 
standard interference term with the weak neutral current the resulting cross section 
for e + p —> e + X in the region of the observed excess at HERA is then determined 
to within about 10 %. 




FIGURE 8. Comparison of the data on the charm structure function with the analysis of ref. 
[10]. The HERA data are from refs [3,4], the Zeus points denoted by stars. The large x data are 
from EMC [17]. 


RE-SUMMED ln(l/X) TERMS 

So far I have discussed fits using LO and NLO resummation of large logQ 2 
terms. Now we should ask is this really good enough when, at HERA, we explore 
x below 1CU 5 , i.e. In 1/x > 9, and the BFKL resummation of In \j x terms i.e. 
Set” (In l/x) m should somehow be included. Catani and Hautmann [24], by incor¬ 
porating kx —factorisation into the collinear factorisation framework, calculated the 
relevant anomalous dimensions to lowest order in a s for each power of ln(l/x) and 
similarly for the coefficient functions. This prompted attempts to calculate struc¬ 
ture functions within the conventional R.G. framework but including the leading 
ln(l/x) terms. Comparisons with data suggested no improvement over fits with 
soft inputs and without ln(l/x) resummation [25,26] and indeed could worsen the 
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FIGURE 9. ^(e.m.) at x = 0.45. The data shown are from BCDMS [14] and from SLAC [15]. 


agreement with data. Furthermore there was a strong dependence on the choice of 
factorisation scheme in which the leading In 1/x terms were included [25-28]. 

The key to incorporating the resummation contributions is to derive an expansion 
for the structure functions which is leading order in both ln(l/x) and a s (Q 2 ) 
and which is renormalisation scheme consistent. Thorne [5] has shown that this 
automatically leads to results which are factorisation scheme invariant and provides 
a framework in which the ‘physical’ anomalous dimensions of Catani [29] emerge 
naturally. This framework is discussed more fully in Thorne’s talk [30]. There is a 
decided improvement in the description of the F 2 data as a result of including the 
ln(l/a;) contributions in this theoretically consistent manner, especially at small x. 
Fig.10 shows a comparison with the low x data and such a fit demonstrating the 
better agreement at very low x. 

Perhaps the sharpest prediction of this approach is the behaviour of F L {x, Q 2 ) at 
small x where the magnitude of the longitudinal structure function is much smaller 
than the ‘conventional’ two loop prediction. Indeed as a final note for this talk 
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FIGURE 10. Comparison of the low x data for F-Fx, Q 2 ) with the leading order renormalisation 
scheme consistent fit from ref [5]. 

I would stress the importance of measuring directly Fj y at HERA as a probe of 
dynamics of small x physics. 
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